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> | . In this paper, the study of the gravitational collapse of cylindrically distributed two perfect fluid system 

has been carried out. It is assumed that the collapsing speeds of the two fluids are very large. We explore 
this condition by using the high-speed approximation scheme. There arise two cases, i.e., bounded and 
vanishing of the ratios of the pressures with densities of two fluids given by c s , d s . It is shown that the 
high-speed approximation scheme breaks down by non-zero pressures pi , pi when c s , d s are bounded below 
by some positive constants. The failure of the high-speed approximation scheme at some particular time of 
the gravitational collapse suggests the uncertainity on the evolution at and after this time. In the bounded 
case, the naked singularity formation seems to be impossible for the cylindrical two perfect fluids. For the 
vanishing case, if a linear equation of state is used, the high-speed collapse does not break down by the 
effects of the pressures and consequently a naked singularity forms. This work provides the generalisation 
6JT), of the results already given by Nakao and Morisawa [1] for the perfect fluid. 
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FT . 1 Introduction 

\o 



One of the central issues in the General theory of Relativity has always been the question of gravitational 
radiation. Historically, the first gravitational radiation problem considered was the radiation from a spinning 
rod [2,3]. One of the important issues of gravitational radiation is the gravitational collapse. General Relativity 
suggests that gravitational collapse of massive objects results in the formation of spacetime singularities in our 
universe [4-6]. The known physical laws break down very near to the spacetime singularity. At this stage, 
quantum theory of gravity is applied to describe the physical phenomena in the neighbourhood of the spacetime 
singularity. One of the important issues "whether the spacetime singularities, (far from the region where the 
gravitational collapse occurs) formed in our universe, can be seen by observer or not" has attracted many people. 

Penrose [7] investigated this problem and suggested a cosmic censorship hypothesis. This hypothesis includes 
two versions, one is called strong version and other is called weak version. According to the strong version, 
a spacetime singularity formed by non-singular initial data is not visible from infinity while according to the 
weak version, a spacetime singularity formed by non-singular initial data is invisible for any observer [8] . The 
singularity claimed by the strong version is named as naked singularity, while the singularity claimed by the 
weak version is known as globally naked singularity. Nakamura et al. [9] conjectured that large spacetime 
curvatures move away to infinity in the form of gravitational radiation in the neighbourhood of globally naked 
singularities (if exist). 

The spherically symmetric systems are considered to be very useful [10] to discuss the problem of naked 
singularity formation. These systems are simple and have well-defined physical significance. We shall take 
non-spherical perturbations to this system [11] to study the problem of gravitational radiation. We consider 
cylindrically symmetric system because there is a degree of gravitational radiation and the spacetime singularity 
formed in this system is naked [12,13]. There is a great literature [14-19] available on gravitational radiation 
through cylindrical gravitational collapse. Some numerical work on the gravitational radiation has also been 
done by different people [20-22]. 
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Recently, Nakao and Morisawa [1] have discussed the gravitational collapse of a cylindrical thick shell com- 
posed of a perfect fluid to the case of non-vanishing pressure. Same authors [23] have also specialized this 
work for the complete dust case. These investigations have provided interesting results about the gravitational 
collapse. This paper has been addressed for the high-speed cylindrical collapse of the two perfect fluids. We 
use high-speed approximation scheme by considering perturbation analysis. For this purpose, the collapsing 
speeds of the two perfect fluids are assumed almost equal to the speed of light. Further, we take the deviation 
of the 4- velocities of the two fluids from null as a small perturbation. We investigate the effects of the pressures 
for the large collapsing velocities. It is verified that our results reduce to the perfect fluid case as obtained by 
Nakao and Morisawa [1]. 

The rest of the paper is outlined as follows. In section 2, we give the basic structure of the spacetime with 
whole-cylinder symmetry [24] . Section 3 is devoted to discuss the high-speed approximation scheme for the two 
perfect fluids. In section 4, we interpret the effects of pressures of the two fluids on the high-speed gravitational 
collapse. Finally, section 5 contains summary and discussion. 

2 Cylindrically Symmetric Two Perfect Fluid System 

The spacetime with the whole-cylinder symmetry is given by the line element [24] 

ds 2 = e 2( 7 -V) ( _^2 + dr 2 } + e 2^ z 2 + e -2V jR 2 # 2 5 ^ ^ 

where 7, ip and R are functions of t and r only. Einstein field equations (EFEs) are given by 

R ab - \Rg ab = «T ab , (a, b = 0, 1, 2, 3), (2. 2) 

where R ao is the Ricci tensor, R is the Ricci scalar, T ao is the energy-momentum tensor and k is the gravitational 
constant. These equations for (2.1) take the following form 

i = (R' 2 - R 2 )~ 1 {RR'(ip 2 + tp' 2 ) — 2RRipi/j' + R'R" - RR' 

- K^{R'T\ + RT r t )}, (2. 3) 
7 = -(R' 2 - R 2 )~ 1 {RR(ip 2 + tp' 2 ) - 2RR'ipijj' + RR" — R'R' 

- K V=g(RT t t + R'T r t )}, (2. 4) 

7 - i'^^-^^'y-gTU, (2.5) 
R - R" = -*V^(T* t +T\), (2.6) 

# + = -^V=9{T t t + T\-T* z + T* 4> ), (2.7) 

where dot and prime represent differentiation w.r.t. t and r respectively. 

There is a literature available [25,26] which shows interesting investigation of the combination of two dif- 
ferent fluids, e.g. Hall [25] presented eight different combinations and investigated their physical as well as 
mathematical properties. Here we have taken one of these combination which is more simpler and has obvious 
interesting properties due to perfect fluid combination. This system of two perfect fluid is also used to obtain 
kinematic self-similar solutions of the Einstein field equations [27]. For the sake of simplicity, we have taken 
only two fluids instead of taking arbitrary. 

We consider the non-interacting combination of two (non-zero) perfect fluids where their (unit, time-like, 
future-pointing) flow vectors u and v are assumed non-parallel. If these fluids have densities pi , P2 and pressures 
Pi, P2, the energy- momentum tensor for two perfect fluid system [25] is given by 

Tab = (Pi + Pl)u a U b + (fi 2 + p 2 )v a V b + (Pi + P2)9ab- (2. 8) 

It is assumed not only that the separate dominant energy conditions hold, but that the extra conditions pj > 
hold and so pj + pj > 0, (j = 1, 2). The Segre type is {1, 1(11)} with corresponding eigenvalues —(pi + p<i) — 
e , (Pi +P2) + e, P1+P2 (repeated), where 

2e = -(pi +p 2 + pi + P2) 2 + {(pi +P2 + Pi + p 2 ) 2 

- A{p 1 +p 1 )(p 2 +p 2 )[l-{u a v a ) 2 ]} 1 2 (2.9) 
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and so e > 0. It is obvious that if energy-momentum tensor represents a perfect fluid with p + p > 0, then it 
uniquely determines the fluid flow u a and the quantities p and p. For the combination of two perfect fluids, the 
physical characteristics of the fluid represented by the flow vector, the pressure, and the density are uniquely 
determined. Since T a b uniquely determines its eigenvalues, p\ + P2 and p± + p2 are determined and p\ + p\ and 
P2 + P2 ar e determined to within an interchange of their values. Consequently, for each choice of the values of 
Pi + pi and p 2 + P2, there is a one-parameter family of values for the quadruple (pi, p 2 , Pi, P2) restricted by 
the requirements that < Pj < pj and pj > 0. 

The components of the 4-velocity u a for the fluid 1 of the whole-cylinder symmetry take the form [1] 

u a =m*(1,-1 + C/,0,0). (2. 10) 

Since u a is time- like, U is positive. Using the normalization property, i.e., u a u a — — 1, u l is written as 

u*= . =■ (2. 11) 

y/U(2-U) 

Similarly, for the fluid 2, the 4-velocity will become 

v a = v*(l,-l + V,0,0), (v a v a = -l) (2.12) 

where v t turns out to be 

p-7+1/' 

W* = ; (2. 13) 

y/V(2-V) 

We define the new variables D\, P\ and £>2, Pi corresponding to the two fluids as 

^(pi+PiK Re^(p 1+Pl ) 
Dl - = VU(2~U) = U(2-U) > (2 ' 14) 

- vn2-v) " v(2-v) ' ( } 

« : - IM? <2 ' 16) 

where g represents the determinant of the metric tensor g a b. The non- vanishing components of the energy- 
momentum tensor T a b turn out to be 

V^gT't = e->-^{-D 1 -D 2 + U(2-U)P 1 +V(2-V)P 2 }, (2.18) 

V=9T r t = e^{D 1 (l-U) + D 2 (l-V)} = -V=gT t r , (2.19) 

sf=gT\ = e->-^{D 1 (l-U) 2 + D 2 (l-V) 2 + U(2-U)P 1 

+ V(2-V)P 2 }, (2. 20) 

V=9T Z Z = e^{U(2 - U)Pi + V(2 - V)P 2 } = ^f~gT* V (2. 21) 

Equation for the conservation of energy- momentum, i.e., T fc a ; a = corresponding to the two fluids can be written 
in the component form as 

d u (Di+D2) = -\(DiU + D 2 V)' + ^{U(2- U)P 1 +V(2 -V)P 2 \ 

+ M {2W-7)-^-7)} 

+ \{U{2 - U)P 1 + V(2 - V)P 2 }(i> - 7 " |), (2- 22) 
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D 1 d u U + D 2 d u V = (1 - U)d u D 1 + (1 - V)d u D 2 

+ \[{U{l-U)D l + V(l-V)D 2 \ 

- {U(2-U)P 1 +V(2-V)P 2 }}' 

- ^-{28 u (yj - 7) - Uty - 7)} 

- ^{2a„(V.-7)-^(^-7)} 

+ |{17(2 - U)P l + V(2 V)P 2 }(i + ^), (2. 23) 

where u = t — r is the retarded time, w = t + r is the advanced time and d u is the partial derivative w.r.t. u. 
The first of these equations is obtained from the i-component and the second is obtained from the r-component. 
The z- and ^-components vanish. 

The energy flux of gravitational radiation can be found by considering the C-energy E and its flux vector. 
The C-energy is defined as [24] 

E=\{1 + e-^iR 2 - i?' 2 )}, (2. 24) 



where E = E(t,r).The corresponding energy flux vector J a is defined as 

V=gJ a = -(E',-E,0,0). (2.25) 

where J a satisfies the conversation law. When we make use of Eqs.(2.3), (2.4) and (2.24) in Eq.(2.25), the 
C-energy flux vector can be written in component form as 



e 



~2 7 



V^gJ* = {RR'{4> 2 + V ) - 2RR^' - /s>/=ff(iZ'T t t + RT r t )}, (2. 26) 

K 

— 2"y 

y^gJ r = {RR{^ 2 + ip' 2 ) - 2RR'ipip' - Kv ^g~(R'T r t - RT\)}. (2. 27) 

The remaining components vanish. 

3 High-Speed Approximation Scheme 

In this section we use linear perturbation analysis w.r.t. a small parameter e. This is called high-speed 
approximation scheme. Since we are considering two perfect fluid system, we take the ingoing null limits of the 
cylindrical two perfect fluids. Using D\, D 2 and Pi, P 2 , the energy- momentum tensor can be written in the 
form 

p 3(^-7) 

Tab = — ^—[Drfakb + D 2 lJ b + e 2 ^{U(2 - U)P 1 + V(2 - V)P 2 }g ab }, (3. 1) 

where 

k a = (1,-1 + £7,0,0), (3. 2) 

l a = (1,-1 + V,0,0). (3. 3) 

In the limits U — > + , V — > 0+, the time-like vectors k a and l a become null vectors. Keeping D\, D 2 , P\, P 2 
finite with these limits, the energy-momentum tensor coincides with the collapsing two null dusts, i.e., 



e 3«>-7) 



<ab 



R 

where 



[D^ah + D 2 lJ b ], (3.4) 



k a -> (1,-1,0,0), (3. 5) 

l a - (1,-1,0,0). (3. 6) 
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Figure 1: The figure shows Morgan's cylindrical null dust solution. The null dust is shown in the shaded region, 
where D(v) > for < v < vw ■ The dashed line on t-axis corresponds to the intermediate singularity at which 
an observer suffers infinite tidal force although any scalar polynomials of the Ricmann tensor do not vanish. 
The dotted line on t-axis is the conical singularity. The short dashed line t = r is the Cauchy horizon. 



It follows that the two perfect fluid system is approximated by a two null dust system for very large collapsing 
velocities, i.e., < U -C 1, 0<V<1. We take the deviation of the four- velocities from a null as a perturbation 
and perform linear perturbation analysis. 

It can be found that for the collapsing two null dusts, the solution is given by 

V = 0, (3. 7) 

7 = 1b{v), (3. 8) 

R = r, (3. 9) 

n{D l + D 2 )e< = (3. 10) 

where -fs is an arbitrary function of v. This solution has been investigated by many people [14,16,17]. We take 
this solution as a background spacetime for the perturbation analysis. It follows from Eqs .(3.4) and (3.9) that 
the energy-momentum tensor diverges at the symmetry axis r = if D\ and D 2 do not vanish simultaneously. 
Using Einstein field equations, the same behaviour can also be shown for the Ricci tensor. This gives rise to a 
naked singularity. It is mentioned here that the same results hold for the perfect fluid case [1] with either D x = 
or £>2 = 0. Some previous studies [16] show that although the Kretschmann scalar vanishes, the tidal forces 
felt by freely falling test particles may become unbounded. This indicates that a non-scalar or intermediate 
singularity is formed there [28,29]. The singularity structure can be shown in the figure 1. 

Now we take large collapsing velocities of the two perfect fluids and use linear perturbation analysis. We 
assume that orders of the variables are such that U = 0(e), V = 0(e) and tp — O(e). The variables 7, R and 
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D\, D 2 can be written as 



e 7 = e 7B (l + <S 7 ), (3. 11) 

R = r(l + 6 R ), (3. 12) 

£>i = D B (l + S Dl ), (3. 13) 

D 2 = D B (l + d D2 ), (3. 14) 



where <5 7 , S R and 5d 1 , Sd 2 are of 0(e) and D B is given by 

D B :=-^^. (3.15) 

It is mentioned here that B stands for background spacctime in the perturbation analysis and "f B (v) is an 
arbitrary function of the advanced time v. We expand the variables U, V, tp, 6 1 , S R , Sd 1 , Sd 2 up to first order 
w.r.t. e. Eqs.(2.3)-(2.7) and (2.22) take the following forms respectively. 

V = 2 K D B e^{5 7 -^+ {6Dl+ 2 SD2) 
{UP 1 + VP 2 ) 



2d v (r5 R )} + (rS R )", (3.16) 



D B 

5, = 2KD B e^{6 7 -^+ { ^±^ 

- - 2d v (rS R )} + (rS R )', (3. 17) 

<5 7 - S^' = -^—(UP 1 + VP 2 ), (3.18) 

rS R {rS R )" = 2ne^{(D B - 2P t )U + (D B - 2P 2 )V}, (3. 19) 

V 1>" - -V' = ~ {(^s - 2Pi)(7 + (.Ds - 2P 2 )V}, (3. 20) 

d u (5 Dl +8 D2 +25 7 -2iP) = -1{(U + V)' -*£L(U + V)} 

+ -^{(upi + vp 2 )- -(up 1 + vp 2 )^-} 

D B av 

1 dDB (U + V). (3.21) 



2D B dv 



Using Eq.(3.21), Eq.(2.23) takes the form 



UP^ + VP? 

d u {(D B - 2P X )U + (D B - 2P 2 )V} = — -. (3. 22) 

r 

The first-order expression for the C-energy w.r.t. e becomes 

E=\[l- e- 2 ^ + 2e- 2 ^ {<5 7 - (rS R )'}} . (3. 23) 

o 

From Eq.(3,15), it is clear that 75 is constant in the region where D B = 0, i.e., in the vacuum region. Then it 
follows from Eqs.(3.16) and (3.17) that <5 7 — (rS R )' is constant in the vacuum region. Thus up to the first order, 
the C-energy is constant in the vacuum region. This implies that, up to the first order in e, the C-energy flux 
vector J a vanishes and up to the second-order in e, it is given in component form as 

V~gJ l = -W- 2 + ^' 2 ), (3.24) 

K 

v^gJ r = -—ipip'. (3. 25) 

It follows from here that the C-energy flux vector corresponds to the massless Klein-Gordon field. 
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4 Effects of Pressure on the High-Speed Collapse 

In this section, we discuss the consequences of the pressures on the high-speed collapse. We take energy densities 
pi , P2 and pressures p\ , p 2 of the two fluids positive and define the following quantities 

c s (u,v) = M, (4. 1) 

V Pi 

d a (u,v) = x f^, (4. 2) 

V P2 

where the subscript s stands for speed. When the dominant energy conditions are satisfied, c s and d s are less 
than or equal to the speed of light [30] . In this paper, we assume that c s and d s are such that 

c s < 1, (4. 3) 

d s < 1. (4. 4) 

It will be shown below that the high-speed approximation does not work for the cases of c s = 1 and d s = 1. 
From Eqs.(2.14), (2.16) and (4.1), it follows that 



ft - j^ft. (4. 5) 



Also, Eqs.(2.15), (2.17) and (4.2) yield 



P 2 = ~^D 2 . (4. 6) 

1 + d s 



Using Eqs.(4.5) and (4.6) in Eq.(3.22), we get 



fl !( l-c s \ TT , f l-d s \ , n Uc s 2 Vd s 2 
l + c s z l + d s (l + c s z )r (l + d a )r 



Integrating this equation, it turns out 

1-c, 2 1-ds 2 f u 2Uc s 2 (x,v) 

[ l + c s 2)U+[ l + d s 2> J 0iv) {l + Cs 2 (x,v)}(v-x) aX 

+ f dx, (4.8) 

Jc7(d) {1 + d s (x,v)}(v - x) 

where U(v) is an arbitrary function of the advanced time v. Eq.(4.8) can be re- written in the following two 
forms either 

,1 + cA.l-dA., .l + c s 2 ..r 2Uc s 2 (x,v) 

- ' -dx 



1 - Cs 1 + d s 1 - c s z J UM {1 



f 

Jui 



u(v) \i + c s 2 (x,v)}(v-x 
2Vd s 2 (x,v) 



u( v ) {1 + d s 2 (x,v)}(v - x] 



dx], (4. 9) 



or 



V l-4 2 l + c s 2 ' 1-4 2 7&(«) {l + c s 2 (x, W )}( V -x) 

Ju(v){l + d s {x,v)}(v-x) 

It can be seen from Eqs.(4.9) and (4.10) that the velocity perturbations U, V are not defined for c s = 1 and 
rf s = 1 respectively. Thus the high-speed approximation scheme ceases to hold for these values of c s , d s . 
Now we consider the following two cases: 

(i) When c s and d s are bounded below by some positive constant; 

(ii) When c s and d s vanish for pi — > oo, P2 — » oo. 
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4.1 Bounded Case 

Here we take c s and d s as bounded below by some positive constants e and / respectively. It follows that the 
following inequalities hold 

2r 2 

" " > e\ (4. 11) 



1 + c, 2 
2d 2 



. > f, (4- 12) 

1 + d s 

where e and / are non-vanishing constants. Using Eqs.(4.11) and (4.12) in Eq.(4.7), we obtain 

1-c 2 1 - d 2 Ue 2 V f 2 

du{{—^- 2 )U + (i— ^ )V] >— + 1*-. (4. 13) 

l + c s 2/ l + d s z v-u v-u 



Integrating this inequality, we have 

(\-^)U + (^-^)V>ln{constxr-( Ue2 + v ^} 
1 H~ c s 1 -\- cb s 

— \n(v - u)du + f — \n(v - u)du, (4. 14) 

where the constant is positive. This shows that the perturbation velocities U, V diverge at the symmetric axis 
r = 0. Consequently, the high-speed approximation scheme breaks down in this case. 

4.2 Vanishing Case 

In this case, we investigate the behavior of the two perfect fluids when c s , d s vanish for the high energy density 
limits, i.e., p\ — > oo, p 2 — > oo. The energy densities will become larger and larger when the fluid elements 
approach the background singularity, i.e., r = 0, (the region of Db ^ 0). Thus the asymptotic behaviors of c s 
and d s turn out to be 

c s 2 (u,v)~C s 2 (t)r« =C s 2 (v-r)r<> ~ C s 2 {v)r\ (4. 15) 

d s 2 (u, v) - D 2 {t)r q = D s 2 (v - r)r q - D s 2 (v)r q , (4. 16) 

where C s , D s are arbitrary functions and q is positive constant. It is mentioned here that D s is different from 
Di, D 2 defined in section 2. Using the above asymptotic forms of c s and d s in Eq.(4.8), we get 

l + C s (v)r q 1 + D s (v)r q q 2 

_ Z_Dl^ [vMrq _ v{0(vlv}{ v_U )q] 



+ -iC ^(x,v)ln{l + C 2 (v)( V —^y}dx 

q Ju(v) ax z 

+ / ( X , v )ln{l + D s 2 (v)( — )"}dx}. (4.17) 

Ju(v) dx 2 



<U(v) 

From here we see that the perturbation velocities U, V remain finite at the background singularity, i.e., r = 
0. Consequently, the high-speed approximation scheme works well in this case. This shows that the naked 
singularity forms at the background singularity through this cylindrical gravitational collapse. 

Now we consider the asymptotic equations of state which have the same asymptotic behavior as Eqs.(4.15) 
and (4.16). From Eqs.(2.14) and (2.15), it follows that the energy densities pi, p2 take the following forms 

2D B V 

p 2 2— 4. 19 
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at the background singularity. Adding Eqs.(4.18) and (4.19) and use Eq.(4.17) at the background singularity, 
we obtain 

or) i 

Pl + 92 ~ ^ + V ) = (function of vf x -. (4. 20) 

When we make use of Eqs.(4.1), (4.2), (4.15), (4.16) in the above equation, it follows that 

Pi + P2 = c s 2 pi + d s 2 p 2 ~ const x ( Pl + p 2 y- q (4. 21) 

where constant is positive. This yields a linear equation of state. Thus the high-speed collapse does not break 
by the effects of pressures if linear equation of state is used. This forms a spacetime naked singularity. 

4.3 Generation of Gravitational Radiation 

It can be seen from Eq.(3.20) that the right hand side represents the source term. This source term contains 
perturbation velocities U, V and if these velocities become larger, also becomes larger. Due to these large 
perturbations of U and V, the C-energy flux in the the asymptotically flat region, given by Eqs.(3.24) and 
(3.25), implies that large amount of gravitational radiations emit. 

The bounded c s , d s case indicates that the perturbation velocities U, V become unbounded at the back- 
ground singularity. For the high-speed approximation scheme, the perturbation velocities are meaningful only 
less than unity. This shows that a significant amount of gravitational radiation is emitted by the pressure 
deceleration at the background singularity. Our result agrees with Piran's numerical results [20] which indicate 
that the large amount of gravitational radiation is emitted when the pressure bounce occurs. 

It has been found that for vanishing c s , d s case that the perturbation velocities U, V remain finite even at the 
background singularity. This result is in agreement with the dust fluid case for the generation of gravitational 
radiation [23]. 

5 Outlook 

Gravitational collapse is one of the most striking phenomena in gravitational physics. The cosmic censorship 
has provided strong motivation for researchers in this field. It is interesting to investigate the gravitational 
waves from a naked singularity formed by the high-speed gravitational collapse by using different symmetries. 

This paper is devoted to investigate the gravitational collapse of the cylindrically symmetric two perfect 
fluids using the high-speed approximation scheme. We have considered two cases, i.e., bounded and vanishing 
of c s , d s . It is shown that the high-speed approximation scheme breaks down by non-zero pressures pi , p 2 when 
c s , d s are bounded below by some positive constants. The failure of the high-speed approximation scheme at 
some particular time of the gravitational collapse suggests the uncertainity on the evolution at and after this 
time. In the bounded case, the naked singularity formation seems to be impossible for the cylindrical two perfect 
fluids. For the vanishing case, it is shown that in the limit of high energy densities, p\ — > oo, p 2 — > oo, the 
high-speed collapse continues until a naked singularity forms. Alternatively, if a linear equation of state is used 
such that pi + P2 is proportional to (pi + P2) 19 with < q < 1, the high-speed collapse does not break down 
by the effects of pressures and consequently a naked singularity forms. 

It follows from Eq.(3.20) that the large amount of gravitational radiation can be measured with the increase 
of the velocity perturbations U and V. This means that the cylindrical two perfect fluids with c s and d s bounded 
below by some positive values will emit more gravitational waves in its collapse due to the unbounded increase 
in U and V which breaks down the high speed approximation scheme. However, for the vanishing of c s , d s 
in the high energy density limit, U, V are finite even at the background singularity. The vanishing of c s , d s 
would generate the larger amount of gravitational radiation keeping the perturbation velocities finite. Also, 
high speed approximation scheme remains valid here. This shows the similarity of the behavior with the dust 
case and consequently generates the larger amount of gravitational radiation. It is mentioned here that if we 
take pi = = pi or p 2 = = p 2 all the results reduce to the cylindrical perfect fluid case [1]. It would be 
interesting to explore this problem by taking some other symmetry. 
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